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Abstract. Constructions of Hadamard matrices from smaller blocks is a well-known tech- 
nique in the theory of real Hadamard matrices: tensoring Hadamard matrices and the clas- 
sical arrays of Williamson, Ito are all procedures involving smaller order building blocks. We 
apply a new block-construction for order 6 to obtain a previously unknown 2-dimensional 
^ family of complex Hadamard matrices. Our results extend the families Dgit) and Bq{9) 

found by various authors recently [T], [1]. As a direct application the existence of a 2- 
^1 parameter family of MUB-triplets of order 6 is shown. 
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ri 1. Introduction 

^ Constructions of complex Hadamard matrices of small orders was originally motivated by 

, G, a question of Enflo, who asked whether for prime orders an enphased and permuted version 

of the Fourier matrix is the only circulant matrix whose column vectors are bi-unimodular. It 
^ was known that this is true for p = 2, 3 but a subsequent general construction due to Bjorck 

[3] (see also the papers written by Munemasa and Watatani [12], and by de la Harpe and 
Jones [6] ) showed that there are inequivalent examples already for any prime p > 7. The only 
^ remaining case p = 5 was settled by Haagerup who has fully classified complex Hadamard 

^ J matrices up to order 5, and showed that Enfio's hypothesis is still true for p = 5. Haagerup 

I also pointed out some possibilities for parametrization in composite dimensions, and in- 

^ troduced an invariant set in order to distinguish inequivalent complex Hadamard matrices 

from each other [5J. Currently the smallest order where full classification is not available 
is order 6. Another significant paper on complex Hadamard matrices was the cathaloge of 
^ complex Hadamard matrices of small orders by Tadej and Zyczkowski who, besides introduc- 

ing another invariant, the defect, listed all known parametric families of complex Hadamard 
matrices up to order 16 [T7]. Most of the presented matrices could be obtained via Diva's 
general method j^, but matrices due to Bjorck [3], Nicoara, Petrescu [13] and Tao [18] have 
also been exhibited. Recently the online version of this cathaloge has significantly been 
extended by new matrices in at least two different ways: firstly a new general construction 
of Butson-type matrices (i.e. matrices built from roots of unities) was discovered by Ma- 
tolcsi, Reffy and Szollosi [10], who used a spectral set construction from [9], while another 
independent construction of Szollosi showed how to introduce parameters to real Hadamard- 
and real conference matrices to obtain parametric families of complex Hadamard matrices 
[TB] . Secondly, new order 6 matrices were constructed by Beauchamp and Nicoara [T] and 
by Matolcsi and Szollosi [TT]. In particular all self-adjoint complex Hadamard matrices have 
been classified, and a family of symmetric matrices has been introduced, respectively. On 
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the one hand, constructing complex Hadamard matrices of order 6 is interesting of its own as 
currently this is the smallest order where full classification of Hadamard matrices is not avail- 
able. While recent numerical evidence suggests that the set of complex Hadamard matrices 
of order 6 forms a 4-dimensional manifold [12] , it seems that describing all of them through 
closed analytic formute remains elusive. On the other hand, complex Hadamard matrices 
play an important role in the theory of operator algebras [H], and also in quantum infor- 
mation theory [20] . In particular, the question whether there exist d + 1 mutually unbiased 
bases (MUBs) in is equivalent to the existence of certain complex Hadamard matrices, 
as such bases can always be taken as a union of the identity operator and a set of (rescaled) 
complex Hadamard matrices whose normalized product is also a (rescaled) Hadamard. For 
a survey on the MUB problem see e.g. [2j, while for a comprehensive list of applications we 
refer the reader to the recent book of Horadam [7j . 

Our paper is organized as follows: in section 2 we derive a two-parameter family of complex 
Hadamard matrices of order 6 by considering circulant block-matrices of order 3. In section 

3 we discuss some connections between this new family and some other previously known 
examples of Hadamard matrices. In particular, we show that besides some well-known 
matrices such as Cq and the members of the generalized Fourier families ^6(1, 3), Fj(l,3), 
the whole affine family DQ{t) and all self-adjoint Hadamard matrices of order 6, denoted by 
Bq{6), belong to our family. All the mentioned matrices can be found online at [IH]. In the 
last section we recall a construction of Zauner [21] to prove the existence of a two-parameter 
family of MUB-triplets of order 6. The main ingredient to his construction is essentially a 
2-circulant complex Hadamard matrix. 



2. The construction 

The main idea of our method is to consider Hadamard matrices with a "highly sym- 
metrical" block structure. Such restrictions made on the matrix implies that "almost all" 
orthogonality conditions immediately hold. We begin our construction with the following 
2p X 2p matrix consisting pxp blocks of matrices A, B and their adjoints A* , B* respectively. 



(1) 



H 



A 

B* 



B 

-A* 



In order to if be a complex Hadamard, one must exhibit certain unimodular matrices A, B 
satisfying the following conditions: 



(2) 



AA* + BB* = 2plp 



(3) 
(4) 



B*B + A* A = 2plp 



AB- BA = 0, 



where /„ is the identity matrix of order p. 



Observe that if we choose A and B to be 
circulant matrices they will commute, and therefore ^ will hold identically, while ^ will 
be equivalent to Hence, by considering p = 3 the building blocks of H can be taken as 



(5) 



A 



a 


b 


c 




' d 


e 


f ' 


c 


a 
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,B = 
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c 
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e 


f 
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and so we have H and its dephased form, Xg 



(6) 
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Xfi 
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ae 


fee 
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cd 
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bd 


be 
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-1 
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bf 


ce 
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-1 
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bd 
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bd 




cf 
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af 
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-1 


be 


cd 


cd 


be 





Let us recall that two complex Hadamard matrices, H and K, are called equivalent, if there 
exists Di, D2 unitary diagonal and P, Q permutational matrices, such that H = D1PKQD2. 

Remark 2.1. Clearly, we are free to set a = c? = 1 by natural equivalence. Also, observe that 
whenever 6 = c we get a self-adjoint matrix by construction. The same holds for the case 
e = / too, as the role of the blocks A, B are symmetric under the equivalence. 

As we have imposed the circularity conditions on the building blocks of if, ([2]) is the only 
equation to be satisfied. In other words, it is necessary and sufficient for if to be a Hadamard 
matrix to find unimodular complex numbers a, 6, c, d, e, /, such that 



(7) 



abode f 
b c a e J d 







holds. At the first glance it seems that we have so much freedom to choose b, c, e, / to 
satisfy ([T]), however, later we will see that there is a really strong connection between these 
seemingly free parameters. Nevertheless we will fully classify this type of matrices obtaining 
a new 2-dimensional family. Let us denote hj (p[x,y] : T x T — * C the following fundamental 
function of ours: 



ip[x,y] ■= x + y + 



1 

xy' 



Now observe, that we have V^[f,^] = f + c~'~a' Hence to satisfy ([t]) one should look for 
certain x,y,u and w G T, such that for some a G ranyj 



(9) 

and 

(10) 



a 



-a 



hold simultaneously. Therefore we should understand the range of ip and characterize the 
set ran(y9 fl ran(— We recall the following well-known 

Fact 1. <^[x, x] = 2x + ^ is a special plane algebraic curve, a three-sided hypocycloid, called 
deltoid. 

The following is also relatively easily seen. 

Fact 2. For any fixed y^ G T, y9[x,?/o] is a sliding line segment with each end on the deltoid 
and tangent to the deltoid. Therefore ^[x,y\ is the union of all such line segments, i.e. the 
whole interior of the deltoid. 
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Figure 1. The intersection of the two deltoids is the fundamental region 
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Let us denote the intersection of the two deltoids above by D := ran(/) fir an (—(/)). It is clear 
that for any a G D one can define a complex Hadamard matrix in the following way: take 
any value of say x, y. Then we have a = l,b = x,c = xy. Similarly, take a] to 

obtain the values of u, v. We have d = l,e = u, f = uv. In particular, we have 



XQ{a) = XQ{x,y,u,v) 
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In the rest of this section we describe an algebraic way of inverting ip, i.e. how we can 
determine x, y and u, v from a given a e D. Considering the equation ip[x, y] = a, we have 



(12) 



X + y + 



1 

xy 



a. 



After conjugating and using the fact that x,y,ET we have 
(13) 



1 1 

- H \- xy 

X y 



a 



Instead of solving the system of equations (12)-(13) we multiply equation (12) by x 7^ 



and (13) by x 7^ and rather consider their sum and difference respectively. In this way 
the variable y vanishes from the difference and we obtain the following cubic equation for x, 
depending on a. 



(14) 



fa{x) := — ax'^ + ax — 1 = 
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It is important to realize that ?/ is a root of ( 14 ) by symmetry as well. Moreover, if x and y 



are distinct roots of (14) then (12) follows. For a G intD the roots of (14) are distinct, and let 
us denote them by /"i, ^2, rs. Our construction guarantees that two of them are unimodular. 
But as rir2r3 = 1 we conclude that the third root is unimodular as well. Hence one can 
choose X as any of ri, r2, ra, and choose y as any other root. We therefore have 6 choices for 
the ordered pair {x,y). 



Finally, let us substitute —a into (14), and denote the roots by qi,q2,Q3- The method to 
determine the values of m, f is completely analogous to what we have presented for x and y. 

For a G intD we therefore have 6 x 6 = 36 choices for the ordered quadruple {x,y,u,v). 
However an easy automatized calculation shows that all of the emerging matrices Xq{x, y, u, v) 
are equivalent to one of the two matrices XQ{ri,r2,qi,q2) or Xj{ri,r2,qi,q2) (note that a 
complex Hadamard matrix is generically not equivalent to its transpose)]^ On the boundary 



of D, however, it is easy to show that the roots of (14) are r, r and ^ and the two families 
Xei^Ti, r2, qi, 52) and Xji^ri, r2, qi, ^2) are equivalent, and hence in this case all choices of the 
quadruple (x, y, u, v) lead to equivalent matrices. 

Finally we note that for every a G D X^i^a) is stable under conjugation, that is XQ{ri,r2, qi, ^2) 
and Xe(ri, r2, qi, ^2) are equivalent. 

By summarizing the contents of this section, we establish the main result of the paper. 

Theorem 2.2. There exist two previously unknown 2-parameter non-ajfine families of com- 
plex Hadamard matrices of order 6, Xe(a;,|/, and XT(x,y,u,v), described by formula 



(11) and its transposed. For a G D the values of x and y are determined as roots of fa in 



(14), while the values of u and v are determined as roots of f^a in (14), respectively. 



Proof. The construction above lead us to a 2-parameter family of complex Hadamard matri- 



ces as follows. For a G D let ri(a), r2(a),r3(a) denote the roots of equation (14), being set as 



continuous functions of a. For a given a G D one can set x = ri(a) and y = r2{a). Similarly, 



substitute —a into (14) and denote the roots as qi{a),q2{a),q3{a), and set u = qi{a) and 



V = q2{a)- Finally, define X^i^a) = Xq{x, y, u, v) as in formula (11 ). We emphasize again that 
easy permutation equivalences show that all choices of the roots ri-^{a),ri^{a),qij^(a),qi^{a) 
lead to matrices equivalent to XQ{ri,r2, qi,q2) or Xji^ri, r2, gi, 52)- 

The main claim of the Theorem is that this family (and its transposed) has not appeared 
in the literature so far. To show this, recall that with the exception of the Fourier families 
FQ{a, b) and Fgla, b), all previously known families of order 6 contain less than two parame- 
ters. Therefore we only need to exhibit one particular matrix from our family X^la) which 
does not belong to the Fourier families. Such a matrix can be obtained by choosing e.g. 
ao = 1 on the boundary of D. It is easy to show that in this case all choices of (x, y, u, v) 
lead to a Hadamard matrix equivalent to Dq, which is not included in the families FQ{a,b) 
and FQ{a,b). Therefore, by continuity, in a small neighborhood U of = 1 the family 
XQ^a) is disjoint from F^i^a^b) and Fl{a,b). Hence, inside this neighbourhood lA only one- 
parameter curves can possibly produce already known complex Hadamard matrices of order 
6, while generically Xfy^a) is indeed new. □ 

This shows that the family XQ{a) is at least locally new, around = 1. We expect that 
more is true: the family Xq{q) intersects the Fourier family only at a = 0. 



^We thank Ingemar Bengtsson for pointing out that XQ{x,y,u,v) and Xq {x,y.,u,v) are generically not 
equivalent. 
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3. Connections to previously discovered families 



In this section we analyze how the obtained new family of complex Hadamard matrices 
Xq is related to the previously discovered ones, such as Bq, Dq, Fq, Mq and 5*6, respectively. 
In particular, we prove that both the Beauchamp-Nicoara family of self-adjoint complex 
Hadamard matrices and Diva's one-parameter affine family is contained in the orbit of XQ{a). 
Thus our construction in some sense unifies and extends some of the previously discovered 
families. 

We shall denote the standard basis of by e^, i = 1, 2, . . . , 6, which should be understood 
as column vectors. Also, for later purposes let us denote by D[a] the discriminant function 
associated to (14), i.e. let ri,r2,r3 be the three roots of fa and define 

D[a] := (n - r2)^(r2 - r-sfira - = \a\^ + 18\a\^ 



(15) 



27. 



Clearly, D[a] G M, and a G D if and only if D[a] < and D[— a] < 0. Note also, that on 
the boundary of D we have D[a] = or D[— a] = 0. 

We begin our investigation with the center of D, i.e. we consider the case a = 0. We have 
the following 



Lemma 3.1. For a = one choice of (x, y, m, v) in formula ( 11 ) leads to a Hadamard matrix 
equivalent to -F6(l,3). 

Proof. Straightforward computation. □ 

Next we classify the "extremal" points of D. It has six points which are farthest from the 
center, and another six which are closest to it. These points will be called "maximal"- and 
"minimal" extremal points of D. 

Lemma 3.2. a) The six maximal extremal points of 3 can be obtained by choosing 

= 1,2,...,6 



(16) = V-9 + 6v^eHf+'=^ 

and lead to matrices equivalent to Cq. 

b) The six minimal extremal points can be obtained by choosing 

(17) ar" = e'''^'^ = l,2,...,6 
and lead to matrices equivalent to Dq. 

Proof. Straightforward computation. □ 

Somewhat surprisingly it turns out that the whole family D^lt) is included in in our family 
Xgi^a). This was actually first found by Zauner [21]. We have the following 

Proposition 3.3 (cf. Ex. 5.7. from [21]). Let D(t) be a complex Hadamard matrix of the 
form 



(1^ 
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-1 
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— i 


-it' 


-1 



where t & T is an indeterminate. Then D{t) has a 2-circulant representation. 
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Diag (1, it, i/t, 1, t, — 1/t) and D2 = 



Proof. Let us define tlie unitary diagonal matrices Di 
Diag (1, i/t, it, 1, 1/t, —t). Then one gets 



(19) 



DiD{t)D2 



1 
it 

i 

t 



1 

t 
1 

it 



it 

i 

t 
1 



1 



-t 



it 

i 

t 



it 



it 

i 

t 



□ 



Corollary 3.4. All members of the Difa-family DQ(t) have a 2-circulant representation. 

Proof. The family D{t) above is trivially permutation equivalent to the Dita-family D^lt^) 
as listed in [17j. □ 

Next we turn our attention to the family of self-adjoint complex Hadamard matrices. 

Lemma 3.5. On the boundary of 3 all emerging matrices are self-adjoint. 

Proof. Let us suppose that a G dB>. Then we have either D[a] = or D[— a] = 0. Suppose 
that D[a] = 0. We have already seen that in this case the roots of (12) are r,r and ^, 
and as we obtain equivalent matrices we are free to set x = r, y = Then, as we have 
b = X, c = xy, the statement follows from Remark 2.1 The case a] = is completely 
analogous. □ 

It turns out that all complex self-adjoint Hadamard matrices of order 6 have a 2-circulant 
representations. 

Proposition 3.6. Let B be a complex Hadamard matrix of the form 



(20) 



B 
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-1 

— X 

_ 1 
y 
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_ 1 
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y 
z 
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-y 
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—xyz 
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xyz 


xyz 
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— X 


-1 



—xyz 

Then B has a 2-circulant representation. 

Proof. Let us define permutational matrices P = [ei, 64, 62, 65, 63, eg] , Q = [e^, 61,63, 64, Cq, 62], 
and the following unitary diagonal matrices Di = Diag (1, 1/ 1/y, —l/{xy^)) 
and D2 = Diag (l, 1/ Xj z^l'^ , 1, —xyz^l'^, y\fz). Here, \fz denotes the principal cubic root 

of z, and z^l"^ is the (slightly abusive) notation of (v^)^- Now we see that D1PBQD2 is 
2-circulant, in particular 



(21) 



D,PBQD2 
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1 
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1 
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—xyz^l'^ 


y\f^ 
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-1 




_^2/3 




-1 
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□ 



As the elegant characterization of Beauchamp and Nicoara [T] shows, all self-adjoint 
Hadamard matrices of order 6 are equivalent to a matrix described by (20). 



Corollary 3.7. All self-adjoint Hadamards of order Q has the 2-circulant representation. 

We close this section with the following remark: matrices Mg and 5*6 are not members of 
the family X^la). It was explicitly stated in |1T], that Mg and Mg are inequivalent, and 
hence a local neighborhood around the one-parametric matrix Mg avoids the family Xg, 
which is stable under conjugation. Clearly, as Sq is isolated, it cannot be a member of a 
continuous family of matrices. 



4. The existence of a two-parameter family of MUB-triplets in 

Recall that a family of mutually unbiased bases (MUBs) {81,82, ■ ■ ■ ,Bk} is a collection 
of orthonormal bases of C" such that |(e, /)| = l/\/n whenever e & Bi and / G Bj for some 
i j. One can assume that Bi is the standard basis, and hence the coordinates of the 
vectors of all the remaining bases have modulus l/i/ra. In particular, the column vectors 
of the remaining bases — up to a constant factor — form complex Hadamard matrices. 
It is well known that at most n + 1 MUBs can be constructed, and this upper bound is 
sharp whenever n is a prime power. On the other hand, when n is composite, not much is 
known about the existence of mutually unbiased bases. For a quick introduction to MUBs 
we refer the reader to [2], [20] • In this section the existence of a two-parameter family of 
MUB-triplets of order 6 is concluded. The method described here was discovered by Zauner 
[2T] . who exhibited a one-parameter family of triplets earlier]^ Interestingly, the heart of 
his construction was the existence of the infinite family of 2-circulant complex Hadamard 
matrices described by formula (19), which he used as a seed matrix for producing MUB 
triplets. We recall his machinery and apply it to the two-parameter matrix Xg(a). First we 
recall a simple, but extremely useful lemma on the representation of 2 x 2 unitaries. 



Lemma 4.1 (cf. Lemma 5.5. from [21] )• Suppose that M is a 2 x 2 unitary matrix with 
entries a, b, c and d. Then there exists u, v,x,y & T, such that 



(22) 



M 



a b 


1 


c d 


~ 2 _ 



U + V 
(u—v) 

X 



y[u-v) 

X 



□ 



Before proceeding we need to introduce some notations. Let T be a 2m x 2m block matrix 
with m X m blocks as the following: 

' A B 
C D 



(23) 



T 



Further let U, V, X, Y are arbitrary unitary diagonal matrices, and let us define the following 
matrices with the aid of the Fourier matrix as 



(24) 



m 

7 
m 



-XFrr 



1 
71 



UF„ 



UYF„ 



VF^, -VYFrr 



The existence of a one-parameter family of MUB-triplets of order 6 was also discovered very recently in 
[8] by a method completely different from Zauner's approach. 
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Note that as is unitary, so are Zi, and hence also 



(25) 



Z^ ^Zi 



(U + V) Fm 

{U - V)) Frr 



F-^{(U-V)Y)F^ 
F-' [X-'Y {U + V)) F„ 



In [21] Zauner characterized 2-circulant unitary matrices in the following way. We quote his 
result with a sketched proof for completeness. 



Proposition 4.2 (cf. Prop. 5.6. from [21]). T is a 2-circulant unitary matrix with blocks 
A, B, C, D if and only if there exist 2m x 2m (rescaled) complex Hadamard matrices Zi, Z^ 
as in formula (24), such that T = Z^^Z2- 



Proof (Sketch). Suppose that Zi,Z2 are given as above. Clearly Z^ ^Z2 is unitary. Also, 
for every diagonal matrix D the matrix F^^DFm is circulant, and hence T is a 2-circulant 



unitary by formula (25). 



For the converse, suppose that T is an arbitrary 2-circulant unitary matrix. Then one can 
write T as 



(26) 



T 



m 









A B 
C ID 











with diagonal matrices A, B, C, D. It follows that T is unitary if and only if the matrices 
(27) Su - 



Ck dk 



are unitary for every 1 < k < m. Now use Lemma 4A to represent Sk with unimodular ele- 
ments Uk, Vk, Xk, Vk, from which one readily defines the unitary diagonal matrices U, V, X, Y, 



and finally Zi and Z2 through formula (24). We conclude by observing that in this setting 
formulas (25) and (26) coincide. □ 



Proposition 4^ describes how to construct a triplet of MUBs from a given 2-circulant 
complex Hadamard matrix T. Clearly, the assumption that T = Z^'^Z2 = Z\Z2 is a complex 
Hadamard matrix implies that {I,Zi,Z2} is a collection of 3 MUBs of order 2m. Note, 
however, that in order to use this construction one needs to begin with a suitable complex 
Hadamard matrix T first, from which the unbiased bases Zi and Z2 can be constructed. 
Clearly, the newly discovered matrix X^la) is a perfect seed matrix for Zauner's construction. 
In summary, we have proved the following 

Theorem 4.3. There exists a two-parameter family of MUB-triplets of order 6 emerging 
from the family X^la) via Zauner's construction described in Proposition 4-2. 



We conclude our paper by the following observation: it is plausible that our new family 
X^la) intersects the Fourier families only at a = 0. If one could exhibit similar families 
Xg'^(a) for all members F{a,b) of the Fourier families, that would provide a rigorous proof 
of the existence of a 4-parameter family of complex Hadamard matrices of order 6. The 
existence of such a family is strongly indicated by the numerical results of [15]. This could 
possibly lead to a full classification of complex Hadamard matrices of order 6. 
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